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We discuss the theory and phenomenology of the interplay between the massless graviton and its massive 
Kaluza-Klein modes in the Randall-Sundrum two-brane model. The equations of motion of the transverse 
traceless degrees of freedom are derived by means of a Green function approach as well as from an effective 
nonlocal action. The second procedure clarifies the extraction of the particle content from the nonlocal 
action and the issue of its diagonalization. The situation discussed is generic for the treatment of two-brane 
models if the on-brane fields are used as the dynamical degrees of freedom. The mixing of the effective 
graviton modes of the localized action can be interpreted as radion-induced gravitational-wave oscillations, 
a classical analogy to meson and neutrino oscillations. We show that these oscillations arising in M-theory- 
motivated braneworld setups could lead to effects detectable by gravitational-wave interferometers. The 
implications of this effect for models with ultra-light gravitons are discussed. 



1 Introduction 

In recent years the old Kaluza-Klein idea of a world with additional dimensions got a strong impact from 
attempts motivated by string theory to resolve the hierarchy problem 1 1 1. Especially fruitful is the idea of 
implementing warped manifolds, for which the metric of the four-dimensional space-time can depend on 
the additional dimensions. The most popular models using warped geometries are the Randall-Sundrum 
(RS) models |2 3 1 in which our observable Universe is treated as a four-dimensional brane embedded into 
a five-dimensional anti-de Sitter (AdS) bulk space. The particular properties of the AdS geometry allow 
to describe the trapping of the gravitational zero-mode on the four-dimensional brane |3|. Displacing the 
second four-dimensional brane in the bulk, one can propose a resolution to the hierarchy problem 1 2 1 . 

However, warped two-brane models have some important additional features. The most interesting ef- 
fect is the existence of light massive graviton modes, which have already attracted a lot of attention in 
the context of higher-spin theories and in a cosmological context [4 5 1. A particularly intriguing devel- 
opment is the construction of semi -realistic braneworld models in which the Kaluza-Klein gravitons 1 (KK 
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1 In the following we will denote also a classical linearized gravitational mode by the term graviton. 
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gravitons) can have an extremely low mass and can be responsible for the mediation of all 1 6 , 7 1 or a large 
part of the gravitational interaction |8 9|. In these setups, the Veltman-van-Dam-Zakharov discontinuity 
(VvDZ discontinuity) 1 10 1 (for a recent analysis of VvDZ discontinuity at the classical level see II II ) is 
either unobservable due to the setup 1 12 1 or cured by non-linear effects 1 13 1, though it is still under dispute 
if these models exhibit a realistic behavior in the limit of strong gravitational fields 1141 . 

An effect so far overlooked is the mixing between graviton modes of different masses which inevitably 
occurs in higher-dimensional spacetime models. This effect is of little relevance in models with small extra 
dimensions and heavy KK modes, because in these models the massive gravitons will rarely be produced 
and therefore are of no astrophysical significance. In models where gravitons can have masses in the sub- 
eV scale like in the RS models, the mixing between the different modes can lead to interesting effects 
in gravitational-wave astrophysics. Gravitational waves (GW's) will then consist of contributions from 
several graviton modes. This will induce a beat in the gravitational wave which could be detectable with 
gravitational-wave interferometers. The parameters of the oscillations will crucially depend on the size 
of the extra dimension — commonly parameterized by the radion-field — , thus suggesting the notion of 
radion-induced graviton oscillations (RIGO's). This effect was investigated in 1 16 1 and will be discussed 
in full detail in the present paper. 

RIGO's are particularly interesting because their observability is not limited to a certain minimal value 
of the extra-dimensional curvature scale. RIGO's could therefore become the first detectable signature 
of M-theory phenomenology. From the very beginning, the affinity of the RS two-brane setup to the M- 
theory compactification of 1 17 1 (see 1 18 1 for reviews of effective five-dimensional M-theory) was noticed 
l2l . Actually, the RS two-brane model can be viewed as a simplified version of this M-theory solution 
which appears five-dimensional in an intermediate energy-range. The RS model simplifies the effective 
supergravity action of 1171 to an Einstein-Hilbert action and neglects the gauge-field part of the M-theory 
effective action. Due to this relationship, the RS model has been used as an approximation for aspects of 
M-theory cosmology where a full-fledged M-theory calculation is not feasible (see 1 19 1 for an example). 

On the other hand, the RS model deviates in major aspects from the requirements of M-theory: First, in 
contrast to the RS model, M-theory restricts the size of the largest extra dimension to typically 10 -31 m < 
I < 10~ 29 m [20 18 1 (see, however, |22|). This rules out the possibility to detect the extra dimension 
by measuring short-distance deviations from Newtonian gravity. Second, in M-theory the requirements of 
having the standard-model group structure and supersymmetry breaking forces one to identify our Uni- 
verse with the positive-tension brane (cf. 1211 ). This prevents one from obtaining the simple geometrical 
solution to the hierarchy problem by locating our Universe on the negative tension brane of the RS model. 2 
Therefore it seems that one has to choose between either the attractive phenomenology of the RS model 
or the appealing theoretical framework of M-theory, although it remains tempting to conjecture that if M- 
theory realizes large extra-dimensions in nature they should exhibit a similar setup as the Randall-Sundrum 
two-brane model. 

Considering RIGO's one may, however, be able to probe several additional orders of smallness of the 
extra-dimensional curvature radius if only the radion acquires suitable values. Although the required values 
of the radion are still far from the parameter range of typical M-theory models, the search for RIGO's thus 
offers a possibility to narrow down the gap between the phenomenologically accessible parameter space of 
large extra dimensions and the predictions of M-theory. 

In our preceding paper 1 26 1 we have studied the four-dimensional effective action in the RS two-brane 
model. This action is a functional of the two induced metrics and radion fields on the two four-dimensional 
branes. Using the results of 1 161 1271 . we shall here extend the study of radion-induced graviton oscillations 
in the two-brane world of 1 16 1 and investigate their phenomenological consequences. 



2 There may remain some loopholes though, see for example 1231 for a construction which would allow us to live on the positive- 
tension brane and have the hierarchy inherited from the negative-tension brane or 24 25 for M-theory constructions which would 
allow the standard-model fields to reside on the negative-tension brane. 
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The structure of this paper is the following. In the next section we will introduce the RS setup and 
describe the corresponding nonlocal braneworld effective action and the equations of motion in terms of 
a nonlocal kernel and Green function, respectively. In Sec. [5] to [5] we will then set up the theoretical 
framework for studying RIGO's by deriving the localized equations of motion from the nonlocal theory. In 
Sec.0we will present an explicit expansion of the Green function in terms of its residues which leads to the 
equations of motion in terms of orthogonal mass eigenmodes of gravitons. In Sec.|4]we will demonstrate 
how these modes can be extracted directly from the action. We find that there are two sets of eigenmodes 
of the action and explain how they are interrelated. To further elucidate this finding, we undertake in 
Sec. 13 a direct diagonalization of the action in a limit in which the two sets of modes coincide. This 
method was first employed in the preliminary study of RIGO's in |27|. The rest of the paper will be 
devoted to phenomenology: In Sec.[6]we discuss the description of RIGO's and their possible observational 
consequences in the RS model. After reminding the reader of the analogous effect of quantum oscillations 
in Sec. 16.21 we first discuss the phenomenology of RIGO's on the positive-tension brane in Sec. 16.31 
This will lead us, in Sec. 16.41 to consider in particular the gravitational waves produced by a network of 
cosmic strings on the negative-tension brane. In Sec. l6.5l follows the discussion of RIGO's on the negative- 
tension brane. Finally, in Sec. the discussion is extended to so-called bi-gravity braneworld models. 
In the concluding Sec. [8] we summarize our results and argue for the occurrence of RIGO's in a general 
class of braneworlds featuring orbifolded compactifications. A detailed presentation of the diagonalization 
procedure for the kinetic and massive parts of the gravitational effective action used in Sec.|5]is given in 
the Appendix. 

2 The RS two-brane model and the braneworld effective action 

In the RS setup |2|, the braneworld effective action is induced by the five-dimensional bulk spacetime 
as follows. We start with the theory having the action of the five-dimensional gravitational field with 
metric G — Gab(x, y), A — (/i, 5), fi = 0, 1, 2, 3, propagating in the bulk spacetime x A = (x, y), x — 
x M , x 5 = y, and matter fields <fr confined to the two branes £±, 

S[G,g,c(>} = S 5 [G} + S4G,g,c(>}, (1) 

where the bulk-space part of the action is given by 

S ^ G ] = Y^G~ 5 J d 5 xG^ 2 ( 5 R(G)-2A 5 ), (2) 

and the brane-part of the action is 

&[G,0,0] = X;jf d i x(L m (cl>,d<P,g)-g 1 / 2 a ± + ^-[K}y (3) 

The branes are marked by the index ± and carry induced metrics g = g^ v (x) and matter-field La- 
grangians L m (4>, d<fi, g). The bulk part of the action Q is characterized by the five-dimensional gravita- 
tional constant G5 and the cosmological constant A5, while the brane parts (0 carry four-dimensional cos- 
mological constants er± . The bulk cosmological constant A5 is negative and, thus, is capable of generating 
an AdS geometry, while the brane cosmological constants play the role of brane tensions and, depending 
on the model, can be of either sign. The Einstein-Hilbert bulk action (|2ji is accompanied by terms in 
containing the jumps of extrinsic curvatures traces [K] associated with both sides of each brane 1 28 1 . Solv- 
ing the five-dimensional Einstein equations with prescribed values of the four-dimensional metric on the 
branes, we have obtained, in |26|, the tree-level effective four-dimensional action. 

In the RS two-brane setup, the fifth dimension has the topology of a circle labeled by the coordinate y, 
—d < V < d, with an orbifold Z2 -identification of points y and —y. The branes are located at antipodal 
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fixed points of the orbifold, y = y±, y + = 0, |y_| = d. They are empty, i.e. L m {(j),d(j> : g M „) = 0, and 
their tensions are opposite in sign and fine-tuned to the values of A5 and G5, 



G 

P' "~ AirGd' 



A 5 = -^, cr+ = -<T_ = 73777. (4) 



Then this model admits a solution to the Einstein equations with an AdS metric in the bulk (7 is its curvature 
radius), 

ds 2 = dy 2 + e- 2M/l r ] ^dx^dx 1/ , (5) 

= y + < \y\ < y_ = d, and with a flat induced metric 77^ on both branes [2 1. With the fine tuning 10}, 
this solution exists for arbitrary brane separation d. 

Now consider small matter sources and metric perturbations Jab(x, y) on the background of this solu- 
tion |3| |29l BD1, 



ds 2 = dy 2 + e- 2 ^/ l V ^dx^ l dx u + j AB {x,y) dx A dx B . (6) 
Then this five-dimensional metric induces on the branes two four-dimensional metrics of the form 

9% {x) = a±V»v + 7%(?), ( 7 ) 
where the scale factors a± = a(y±) are given by 

a+ = 1, a- = e- d / 1 = a, (8) 

and 7av ( x ) are tr, e perturbations by which the brane metrics g£ v (x) differ from the (conformally) flat met- 
rics a^rjuu of the RS solution (jSJi- The variable a(y) represents the 'modulus' — the global homogeneous 
part of the radion field determining the interbrane separation. 

Due to the metric perturbations, the branes no longer stay at fixed values of the fifth coordinate y. Up 
to four-dimensional diffeomorphisms, their embedding variables consist of two four-dimensional scalar 
fields, the linearized radions ip^(x), and the braneworld action can depend on these scalars. Their geomet- 
rically invariant meaning is revealed in a special coordinate system where the bulk-metric perturbations 
1ab( x , v) satisfy the "Randall-Sundrum gauge conditions" 7^5 = 0, 7 All / l/ = Zi't = 0. In this coordinate 
system, the brane embeddings are defined by the equations 

£± : v = y± + -tV^Oe), y+ = o, y-=d. (9) 



Therefore the radions describe the bending of the branes in the bulk (cf. 1291 \30l ). 

In 1 26 1 we have derived the effective braneworld action in terms of the four-dimensional on-brane 
metrics and nonlocal matrix-valued form factors. We call the part of the action quadratic in h^ u (x) 
— the transverse-traceless parts of the full metric perturbations 7^,(2;) on the branes — the graviton 
sector. It can be written as the following 2x2 quadratic form in terms of the metric perturbations and the 
special nonlocal operator F(d), 



s 1-^— /d 4 zV h-^ {u) 

<->grav [ Huv \ — 1 / " X r, [ a iiv a L 



JgravL ^ J ~ IGTrGj 2 1 ^ ^ P 



(10) 



The effective gravitational constant G4 is given by G4 = G5/I, As was shown in Ref. |26|, the operator 
F(n) is a complicated non-linear function of the D'Alembert operator □ expressed by means of Bessel 
and Neumann functions of the arguments ly/n and ly/n/a, 



F(D) = ' 



ly/UU + (ly/u/a) 



1 ^MVO) 

TT a 



(ID 
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where 



u ± (z) = Y 1 ± J 2 (z)-J±Y 2 (z), 

J+ = ,h{l^u), 
Jr = JiQVa/a), 



Ff = YxilVn/a). 



(12) 
(13) 
(14) 



The action dlQ> should be amended by the standard coupling of the transverse-traceless gravitational modes 
to the transverse-traceless part of the stress-energy tensors, T~,, on the two branes, 



d A x\h + h~ 1 



T+ 



The inverse of the kernel F(n) is the Green function G(n) of the problem, 
G = F _1 . 

With the abbreviations given in G3 - (O, it reads 



(15) 



(16) 



G(D) 



1 



l 2 n J+Yf - Jf Yj 4 " 



2 

7T 



lVou+(ly/n/a) 



(17) 



It has first been calculated in |50|. With the help of G(d), one finds the equations of motion for the 
transverse-traceless sector 



h 



//.;/ 



-8ttG 4 1 2 G(a) 



T+ 
T~ 



(18) 



corresponding to the variation of the combined action of dlOt and d!5l >. In this paper, we study in detail the 
properties of the model in the low-energy limit when ly/n -C 1 but when ly/u/a can take arbitrary values 
in view of the smallness of the parameter a — e~ d l l (large interbrane distance). 

As we will only be considering the transverse-traceless part of the gravitational dynamics, tensor indices 
will be omitted for the rest of this paper. 



3 The eigenmode expansion of the Green function 

Both the Green function G(n) and the kernel of the action F(n) are highly nonlocal. Nevertheless, it is 
possible to obtain a conventional interpretation of the equations of motion and the action in terms of an 
infinite tower of orthogonal Kaluza-Klein modes. We will first discuss the recovery of the KK tower from 
the Green function because this procedure works in a more direct way than the recovery of the particle 
spectrum from the action, which will be discussed in the Sec.|4] The elements G x x (x = ±) of G(n) are 
meromorphic functions of □. If we consider a concentric circle C n around n = so that C n includes the 
first n poles and let its radius R n — > oo as n — * oo, we find the falloff property 

|G XX (D)| < ei?„, i?„->oo, (19) 

for any small constant e and all □ on C n . Therefore we can apply the Mittag-Leffler expansion, which 
provides an expansion of a meromorphic function in terms of its poles, to the Green function in order to 
obtain a representation of G(n) as a sum of scalar propagators. 
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In the limit ly/n <C 1 but ly/o/a arbitrary, the Green function G(n) is given by 3 



G(D) 



2 

Wn 



2 J^ly/o/a] 

1 



111 



ly/O 



W Mly/U/a] 



^J^lVa/a] 
a Miyn/a] 



(20) 



For convenience we use □ as the fundamental variable of the expansion and not ^f^. The Mittag-Leffler 
expansion for an element G x x of the Green-function matrix reads 



Gxx{0) 



□Res[G > 

i=0 



mf)} 



mj[n-m 2 ] 

where the first pole is at □ = 0, i. e. at mo = 0. One can write 



Gxx(n) = ^ 



Res[G xx (n = m?)] 



□ — mf 



E 



Res[G xx (n = 



(21) 



(22) 



As the second sum contributes only a constant, we can drop this part and obtain our final result for the 
elements of the Green function, 



Gxx(n) = ^ 



Res[G xx (n = mf)} 



i=0 



□ - mf 



(23) 



In this way, each nonlocal element of the nonlocal Green function G(d) can be represented as an infinite 
sum of scalar propagators with different masses. The mass-squares are given by the positions of the poles 
of G(n). From the representation of the Green function J20l >. we immediately infer that all elements of 
G(d) have poles at the same values of □. This is a direct consequence of the common prefactor displayed 
in (1171 which contains all poles of G(n). One pole is located at □ = 0, which is also an exact pole of (1 171 . 
An infinite sequence of poles is located at the roots of the Bessel function J\, Denoting the argument at 
the ith zero of J\ by ji, we find that these poles are located at 



□ = (j ia /l) 2 = m 2 , 
where the first few values of ji are given by 

ji« 3.832, j 2 « 7.016, j 3 w 10.173. 



(24) 



(25) 



Whereas the first term of the sum d23t from the pole at □ = will provide us with the effective four- 
dimensional massless graviton, the infinite series of poles at is responsible for the generation of the KK 
tower of massive gravitons. The matrix of residues of G at the pole □ = is given by 



Res[G(n = 0)] 



1 



l 2 (l-a 2 ) 

and the residues for the infinite sequence of poles at □ = mf are found to be 

2a 2 



Res[G(n = mf)} 



I 2 



l/(J2[lm t /a}) 2 -1/Japmi/o] 
— l/Ja[/mi/o] 1 



(26) 



(27) 



3 Note that the falloff property 1191 is not valid for the element G++ in the approximation 1201 but only for the full expression 
1171 . This does of course not prevent us to use the approximation 1201 for an approximate determination of the residues and thus for 
an approximate determination of the Mittag-Leffler expansion. 
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The residues can be factorized as 
Res[G(n = m,)] 

with 

V2 



T 

ViVj . 



v 



1 



V2o 



l/J 2 [lmi/a] 
-1 



i > 1. 



(28) 



(29) 



This property reflects the fact that the residues of the poles of the Green function are projectors on the 
corresponding propagating modes of the theory. These energy eigenmodes will be recovered from the 
action in the next section. First, however, we should display explicitly our newly acquired representation 
of the Green function: 



G(D) 



i 



oo 



2a 2 



^ l 2 {n - m?) 



l/(J 2 [lmi/a)) 2 
-l/(J 2 [lmi/a]) 



-l/(J 2 [lm l /a}) 
1 



(30) 



4 The spectrum and the eigenmodes of the effective action 

The eigenmodes of the equations of motion can be recovered from the action i ll Ob and its kernel F(n) il U . 
We will find that the situation has many similarities with resonance theory, and our discussion will parallel 
many of the considerations of 1 3 1 1 from which we will also adapt our nomenclature of eigenstates. 

In the low-energy limit, ly/n <€. 1 but l^/u/a arbitrary, the kernel of the action F(n) can be approxi- 
mated by 

l/J 2 [lVa/a\ 



F(n) « 



1 



l/J 2 [ly/o/a] — 



2 Ji[U/n/a] 



(31) 



ly/na J 2 [l\/u/a] 

The typical way to extract the particle content from an action with a matrix-valued kernel is its diago- 
nalization in terms of normal modes. However, as found by the application of the Mittag-Leffler expansion 
of the Green function in the last section, the number of propagating modes enormously exceeds the number 
of entries in the 2 x 2-matrix F(n) and thus the modes do not diagonalize the quadratic action fllOi in the 
usual sense. This behavior could also have been anticipated from the nonlocality of F(n). The propagating 
modes are the zero modes of F(d) which solve the matrix-valued nonlocal equation 



F(n)hi(i) = 0. 



(32) 



These modes are conveniently split into a scalar part hi(x), depending on the four-dimensional spacetime 
coordinates, and a spacetime-independent (isotopic) vector part Vj, 



hi(x) = hi(x)vi . 

The eigenvector equation d32i is accompanied by its consistency condition, 

detF(n) = 0, 



(33) 



(34) 



which corresponds to picking the poles of the Green function G(d). The condition i34i yields the mass 



spectrum of the theory given by the roots of this equation, i. e. □ 



Hence the hi(x) above are Klein- 



Gordon modes and the isotopic vectors of the propagating modes v, are zero eigenvectors of F(to 2 ), 

(n-mf)hi{x) =0, (35) 
F(m?) Vi = 0. (36) 
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Thus we obtain the Kaluza-Klein spectrum which contains the massless mode i = 0, too = 0, and a tower 
of massive modes. In the low-energy approximation of F(rj), see Eq. (13 1 1 , its masses = aji/l are 
given by the roots of the first-order Bessel function, J\ (ji) = 0. The isotopic structure of their Vj is given 
by the vectors d29t which were found to factorize the residues of the Green function. 4 In view of standard 
arguments of gauge invariance, the massless graviton ho has two dynamical degrees of freedom, while the 
massive tensor field Km posesses all five polarizations of a generic transverse-traceless tensor field. 

The action i ll Ob is not, however, diagonalizable in the basis of these states because under the decompo- 
sition h(x) = X)j=o hi(x)*Vi (with off-shell coefficients hi(x)) the cross terms intertwining different i's 
are non-vanishing, 



vf F(n)vj + 0. 



(37) 



A crucial observation is, however, that the diagonal and non-diagonal terms of this expansion are linear 
and bilinear, respectively, in on-shell operators □ — mf , 



vrF(n) Vl 



0-i 



dn 



-v, 



,2 



vf F(D)v, = M i:j (D)(a - m?)(n - mj), i £ j, 



(38) 
(39) 



where higher powers of ( □ — mf ) have been dropped in J38i . and My ( □ ) is non-vanishing at both n = mf 
and □ = m|. Therefore, the non-diagonal terms of the action do not contribute to the residues of the Green 
function G(d) of F(n). The normalization of v^, as given by i29\ . automatically yields in J38I a unit 
coefficient of (□ — mf). 

Thus we have found the particle spectrum of the action in terms of modes which turn out to be non- 
orthogonal off-shell but become orthogonal if one considers the action on the mass-shell. In the theory of 
atomic and nuclear resonances, such non-orthogonal energy eigenstates are called Siegert states after their 
introduction in |32| (cf. [31 1). 

This particle spectrum can also be recovered from the two orthogonal eigenstates of F(n). These states 
which will be energy dependent, i. e. not mass eigenstates, have also been extensively used in resonance 
theory and are named Kapur-Peierls states according to their first use in 1331 1341 (cf. also 13 

In resonance theory, Kapur-Peierls states bj(n) diagonalize the matrix Hamiltonian of the system (the 
analogue of our F(n)). For our 2x2 matrix-valued kernel F(n), this property reads 



F(n)b a (n) = A s (n)b s (n), a = 1,2, 
bj(n)b s /(d) = 5 SS ,, 



(40) 
(41) 



where the second equation describes the orthonormality of the Kapur-Peierls states. 

In the limit l\/n <C 1, the energy-dependent eigenvalues of F(d) [as determined by Eq. ( II 1H are found 
to be 



Ai D = — + O(a), A 2 □ = V 1 *' +0 (o . 
2 a J2[imi/a\ 

The corresponding eigenvectors (normalized to unity to leading order in a) are 



(42) 



bi(n) = 



1 

aly/n 
2 Ji [lmi/a] 



b 2 (D) = 



1 



4 The vector vo is here actually given by vo = v^2A[l o?] T , 
because here the Vj's are obtained for the approximation 1311 of F. 



(43) 



i. e. the leading order in a of the exact vo as given by 1291 . 
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where we have introduced 

Z 





2 1" 


a 





Ji[lm,i/a] — J2[lnii/a] 



(44) 



The eigenvalue Ai, given in ( I42t . has only one root at the value □ = 0, whereas the eigenvalue A2 has 
infinitely many zeros at the zeros of the Bessel function J\, i. e. at □ = mf. Approximating A2(d) around 
its zeros by the first term of its Taylor expansion, we find for the eigenvalues 



l 2 n I 2 
AiC) = — , A 2 (n) = ^(n 

At □ = the eigenvector bi takes the value 



mf). 



MO) 



1 



and at energies □ = mf the eigenvector b2 becomes 



b 2 (mf) 



1 

-J 2 [lmi/a] 



(45) 



(46) 



(47) 



As expected, the mass levels of the Siegert states equal the zeros of the Kapur-Peierls eigenvalues. The 
massless graviton Siegert state corresponds to the single root of the eigenvalue Ai(d), and the tower of 
massive Siegert states arises from the infinitely many zeros of the eigenvalue A2(d). It is easy to find 
the relation between the Siegert eigenvectors and the Kapur-Peierls eigenvectors by taking from Eq. J40i 
that the vectors b s (n) are eigenvectors also of F _1 (n) = G(n) with eigenvalues 1/A s (n). Using the 
completeness and orthonormality of the b s (n)'s, we obtain a representation of G(n), 



G(D) 



E 

8=1,2 



b,(D)tff(D) 
A S (D) 



As the poles of the Green function lie at X s (□) =0 we find for the residues of G(n) 

b s (n)bj(n) 



lim (□-mf)G(n) 



d\ s (a)/dD 



(48) 



(49) 



where s = 1 for mo 
we obtain 



and s = 2 for m., with i > 1. Comparing this expression for the residues with 



v = 



dXi(p) 
du 



-1/2 



bi(n; 



□=0 



( dXxju) 
\ da 



-1/2 



ba(n; 



i > 1. 



(50) 



The decomposition of the action into Siegert eigenmodes (I33t provides us with the conventional particle 
interpretation of the propagating modes of the nonlocal operator F(n). It clarifies their role in its Green 
function ( I30l > that mediates the gravitational effect of matter sources. Amended by the matter action on the 
branes, the effective action of the graviton sector reads 



S[h%} = 



d 4 x 



32ttG 4 I 2 2 



(51) 



where T is the column vector of the transverse-traceless part the stress-energy tensors on the branes and h 
is now given by the sum of Siegert modes (1331 . 



(52) 



i=0 
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Varying this action with respect to each Siegert mode and decomposing the results to recover the equa- 
tions of motion for h, we obtain the linearized equations of motion. Their solution 



87rG 4 / 2 G ret T (53) 




is expressed in terms of the retarded version of the Green function J30i . We note the explicit expressions 
for h + and h~ obtainable from J50i and (1521 : 

h (x) + a y2 xr/ 1 / i ^) . (54) 
^ J2[lm i /a\ J 

hp(x) — aj^hjfa) J . (55) 

H / 

It is particularly interesting to study the asymptotic behavior of the expressions J54t and J55I . In the 
long-distance limit a — > 0, we have 

h+ ~ ho + 0{a) , (56) 
h~ a 2 ho — a J^ ;>1 hi . (57) 

Thus, in this limit, the /i + -mode practically coincides with the massless graviton ho, and the mixing with 
the other modes is suppressed by the factor a. On the other hand, the S_-brane mode h~ is almost 
exclusively composed from the massive modes, and mixing with the massless mode is suppressed, whereas 
the mixing between the massive modes can be called maximal, i. e. h~ is proportional to the sum of the 
massive modes. Moreover, the contribution of both the modes hi and ho to the /i~-mode in this regime is 
suppressed by the common factor a. There is also a regime of maximal mixing between all modes for the 
metric perturbation h~ on the £_-brane at a = 1/ \[2 where 

h~ ~ h a + V hi. (58) 
^ — 'i>l 



5 The graviton effective action in the diagonalization approximation 

In the last section we have found that the Siegert state representing the massless graviton corresponds to 
the zero of the eigenvalue Ai of the kernel of the action F(d), whereas the Siegert states corresponding 
to all of the KK modes stem from the zeros of the second eigenvalue A2. This opens the possibility of a 
low-energy limit in which the Siegert and Kapur-Peierls states coincide. In the approximation 

\/D<l, l/aVn<A, (59) 

one can use the small-argument expansion for all Bessel and Neumann functions appearing in F(n). If we 
truncate this expansion at £?(□), we obtain an expression for the kernel F(d) which is entirely local. The 
choice of the second part of the low-energy condition (I59> guarantees that the approximation is still valid 
slightly above the mass of the first KK mode, 

mi = jia/l (mi/asi 3.83 (60) 

(cf.!25>. In fact, at least for Bessel functions of the second kind, the small-argument expansion is considered 
as the method of choice in numerical calculations for arguments < 4, see [..351 ■ However, we should be 
prepared for some deviation of the expansion-result for the first KK mode from that of the previous sections 
due to the low order of truncation of the expansion. Taking this into account, the result will turn out to be 
surprisingly accurate. 
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In the low-energy approximation ( I59l l. the kernel F/Z 2 can be represented in the form |26 



F(D) 
Z 2 



= -M + Dn +o(d 2 ), 



m =pT 



D 



,,2 



3(1 + a 2 ) 5 



-a 

a 2 + 3 
2 



-a 
1 



(61) 
(62) 

(63) 



The operator D represents the kinetic part of the action (llOi . while the operator M plays the role of the 
mass term. Calculating the trace and the determinant of the matrix D, it is easy to see that it is positive 
definite for any value of the parameter a. The determinant of the matrix M is equal to zero, therefore M is 
degenerate. Using the positive-definiteness of the matrix D, one can diagonalize both the matrices D and 
M simultaneously. 

Our purpose here is not simply to diagonalize these matrices, but also to make a canonical normalization 
of the kinetic terms. With this choice, the action for gravity eigenmodes with fixed masses, representing 
some mixtures of the functions h + and h~, will have the form 



<5graviton[^0j /Im] 



16 Iq 4 J d * x (jhoDho + j/ijw(d - M 2 )h M ^j , 



(64) 



where M is the mass of the massive gravitational mode. 

The details of the procedure of diagonalization are presented in the Appendix. Here we shall give the 
final results. The expression for the mass of the massive graviton mode Km is 



Af2= 24a 2 (l + a 2 ) 



(65) 



Z 2 (l-a 2 ) 2 ' 

At small values of the parameter a, i. e. at large distances between the branes, the squared mass behaves as 



, 24a 2 

M 2 ~ — (66) 

and tends to zero together with a 2 . This result is to be compared with the more accurate expression for the 
first KK mass in relation J24I . The large discrepancy is easily explained from the fact that taking the limit 
a — > is strictly speaking not legitimate in our approximation because the second condition of ( 15 9t would 
not be fulfilled anymore. 

The relations describing the transition from the old graviton variables h + and h~ to the modes hu and 
ho have the following form, 

h M = -=- [ah+ , (67) 

V3(l + a 2 ) V a J 

ho = ^-^(h + + h-). (68) 
1 + a z 

The inverse transformations are 

h+ = (ho + V3ah M ), (69) 



VI- a 



h~ = -j=^= (ah Q ~V3h M ). (70) 
VI — a 
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The corresponding expressions for vo and vm are 



Vo 



V2 



1 



l\/T 



1 

-1 



(71) 



which should be compared with the expressions for vo and vi in d29i . As might have been anticipated 
from the fact that the residues depend only on the (— l)th term of the Laurent expansion of the function 
considered, we find exact agreement for Vo. 

Amending the action (1641 by the matter action \\5\ . we can obtain the equations of motion for the 
metric perturbations h + and h~ on the brane by varying the total action with respect to ho and Km and 
then combining their equations of motion according to the transformations ( I69> and ( 17 Oi l. 



h+ 
h~ 



= -16ttG 4 



1 


1 


'la 2 " 


□ 


1 - a 2 

ret 1 " 


a 2 a 4 _ 



1 



□ - M 2 



3a 2 



1 



1 

-1 



1 



T+ 



(72) 



In the language of the Siegert and Kapur-Peierls eigenmodes of Sec. 0] we have chosen for the diago- 
nalization an approximation in which the Siegert eigenmodes become orthogonal and the Kapur-Peierls 
eigenvalues and eigenstates energy-independent. Thus, in the approximation ( 15 9> . the Siegert and Kapur- 
Peierls descriptions coincide. The first Kapur-Peierls state coincides with the Siegert state of the massless 
graviton, while the massive Siegert states become a kind of "collective" massive state which coincides with 
the second Kapur-Peierls mode. 



6 Graviton oscillations 



6. 1 Gravitational waves on the S + -brane 

After the extensive theoretical considerations of Secs.|4]and|5] we will now return to the equations of mo- 
tion J 1 81 and study their phenomenological properties. We will proceed to show how graviton oscillations 
arise naturally in the propagations of gravitational waves in this model. For studying the properties of 
graviton oscillations, it is convenient to abandon the condensed matrix notation of Eq. dl 81 and to study 



the behavior of metric perturbations on the S + -brane and £_ 
Using the spectral representation ( I30> . we find for h + , 



-brane separately, starting with the £ + -brane. 



16ttG 4 



□ 



(T+ + a 2 T-) 



OO 

E 

i=l 



16ttG 4 



□ 



Ji 



T+ 



-T~ 



(73) 



where we have introduced J<z = J2[lmi/a] ~ 0.403. We consider astrophysical sources of equal intensity 
at x = on both branes with a harmonic time dependence 



(74) 



where /i is the mass quadrupole moment of the source. The sources with an oscillating quadrupole moment 
will generate gravitational waves. At a distance r from the source, the waves on each brane are given by a 
mixture of massless and massive spherical waves. On the S + -brane, this superposition is given by the sum 
of the contributions from the sources on the E+- and £_-brane, respectively, 



h+[T+] = Ae-*"* ( e^ r + V j-rr^—^ 

h + \T-} = Aa 2 e~ wt ( e l " r - V 1 V" 2 -"^ | 
^ f-lMlmi/a] J 



(75) 
(76) 
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where A = 4 G4 fi/r is the amplitude of the massless mode generated on the E + -brane. If the frequency uj 
with which the source is oscillating is smaller than the mass nii of a KK mode, u> < rrii, then the exponent 
of the corresponding mode will be real and negative, instead of imaginary, and the mode will be decaying 
instead of oscillating. In effect all modes with masses above the frequency of the source will have died off 
after propagating for a short distance and only those modes with a mass below the frequency of the source 
will contribute to the gravitational wave in the far field. Thus, as long as we are interested in GW's, we can 
content ourselves with including only the terms with mi < uj into the sums of Eqs. (175 \ and i76\ (cf. l6l). 

In order to have a simple situation, we consider a source with a frequency above the mass threshold 
of the first massive mode but below the threshold of the second KK mode, i.e. mi < u> < m?,. Then 
only the massless and the first massive mode are excited and produce long-range gravitational waves. At 
a distance r from the source, the waves on each brane are given by a mixture of massless and massive 
spherical waves. On the S + -brane this superposition is given by the sum of the contributions 



The above results have interesting physical consequences. We have seen that the observable transverse- 
traceless metric on each brane is indeed a linear combination of different gravitational modes with fixed 
masses. One of these modes appears to be massless, while the others are massive with their masses de- 
pending on the distance between the branes. 

This situation is quite analogous to the well-known mixing arising in the context of kaons 1361 1371 or 
other bosons [38 1 or with neutrino oscillations |39|. As a matter of fact, the phenomenon which we en- 
counter here is not quantum oscillations but classical ones. Nevertheless, the interference effects between 
the different gravitational modes, arising in the process of their propagation in space (on the brane), are 
analogous to those of quantum oscillations. This is quite natural because the oscillations of quantum parti- 
cles are totally conditioned by their wave nature (see, for example, Ref. |37|), and both effects are caused 
by the superposition of modes with different dispersion relations. 

However, there is one important difference: it does not matter for quantum oscillations if one considers 
the superposition to consist of two components of different energies but with the same momentum (as 
is usually assumed) or if the components are supposed to have the same energy but different momenta 
1401 . For a mixture of classical gravitational waves, however, it is natural to assume that they have been 
generated by the same source oscillating with a frequency to. Therefore one should consider oscillations 
between modes of the same frequency to which have different wave-numbers k, and not vice versa. 

Let us recall briefly some basic formulas from the theory of quantum oscillations. Suppose that a 
certain particle can exist in two different observable states \ipi) and ^2), which correspond to a quantity 
conserved in some interaction but which do not have fixed masses. For example, kaons can be created or 
detected in states with fixed strangeness, which are not eigenstates of their Hamiltonian, and according to 
the theory of neutrino oscillations, which is now supported by good experimental evidence |41 1, the states 
of the electron, muon and tau neutrinos are also mixtures of different Hamiltonian (i. e. mass) eigenstates. 
Generally, this situation can be described for the two-state case by 




(77) 



(78) 



from sources on the £+- and £_-brane, respectively. 



6.2 Quantum oscillations — an analogy 



\tpi) = cos0|V>m) +sbx9\ip m ), 
\ip 2 ) = - sm9\ip M ) +cos8\ip m ), 



(79) 
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where the states \i/jm) and \ip m ) are states with masses M and m, respectively, and the angle 9 parameter- 
izes the mixing between the states. If at the initial position of an evolution at x = somebody observes, 
say, a particle in a state 

ma)) = IV'i), (80) 

propagating in space with the energy w, then at the position x the state J80i will become equal to 

=C086e- ikMX \iP M ) +sm9e- lk ™ x \i> m } 

= cos0e-V fc 3+ M2a! |^ M ) +sin0e-V fe §+ m2a: |^ m ) ) ( 8 i) 

where fco = w - Now we are in a position to calculate the probability of observing the particle in the state 
at the position x which is given by the formula 

Pr{x) = |(^|V(*))| 2 = 1 - sin 2 (20) sin 2 ^ . (82) 

Correspondingly, the probability of finding the particle in the state \1jj2) is 

Pz(x) = \m^{x))\ 2 = sin 2 (20) sin 2 ~*™) x . (83) 

It is easy to see that the oscillations of the probabilities P\ and P 2 are strongest when sin (26) = 1, i. e. 
when 

9 = ±|. (84) 

This situation is called the case of maximal mixing between oscillating particles. Looking at the expres- 
sions for the probabilities d82l > and 183t . one can see that they represent periodic functions with a period 



v/fc 2 + AP - v'fcg + m 2 ' 

This is the oscillation length of a quantum state with the energy u> = fco. 

In the case of maximal mixing (I84L detecting the particle at the distance L/2, one finds with 
probability P^L/i) = 1 the particle in the state |^a)- in other words, the particle in the state 
is totally transformed into the particle in the state ^2)- The time-dependence of the oscillation can be 
inferred from this expression by considering semi-classically the propagation speed of the corresponding 
particles. 

The description of quantum oscillations presented above coincides with describing classical oscillations 
of superpositions of waves with different dispersion relations. In this case, instead of the probability of 
detecting a particle of a certain type, one can calculate the amplitude of the wave. 

6.3 Gravitational-wave oscillations on the E+-brane 

For the superpositions d75i and (I76> . the amplitudes detected by a gravitational-wave interferometer are 
given by the absolute values of the waves, 

\h+[T+}\ =A+ 
\h + [T-}\ =A~ 



1 - 



4a 2 J? 



4J 2 



2 • 2 
sin 



/7ITN 



-,1/2 



(J2 1) 



• 2 (Kr\ 

r? sm \T) 



1/2 



(86) 
(87) 
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Here, L is the oscillation length of the amplitude modulation of the gravitational wave (GW), 
2tt 



lu — y 'uj 2 — m\ 



(88) 

This corresponds to Eq. (I85> with one mode being massless. The pre-factors of the amplitudes d86i and 
d87i are given by 



A + =[l + jt)AttA, (89) 
A~ = ( — - 1 ) a 2 A w 1.5 a 2 A, (90) 



where the approximations are valid in the limit a <C 1. We find oscillations in the amplitudes of the 
waves from both sources. For a GW generated by T + , see Eq. d86i . the oscillating part of the amplitude is 
suppressed by a factor of a 2 compared to the constant part in the limit of large brane separation, a <C 1. 
The amplitude of the GW produced by T~, see Eq. (I87> . is totally oscillating, regardless of the inter-brane 
distance. 

In the limit of a source-frequencyjust above the mass-threshold co > mi, the oscillation length (18 8 1 of 
these radion-induced gravitational ("graviton") wave oscillations (RIGO's) tends to 

£>-• (9D 

This limit is mainly useful to set a lower bound on the oscillation lengths as in this limit the propagation 
speed of the massive mode will tend to zero. Therefore, in this limit, the massless and massive wavetrains 
from a source would soon become spatially separated from each other and, far from the source, one would 
detect two separated wavetrains instead of one wavetrain exhibiting RIGO's. The astrophysically more 
interesting limit is the case u> 3> m\ in which the expression for the oscillation length becomes 

L = 4tt^- . (92) 

We can express the minimal oscillation length (19 11 of RIGO's in terms of the AdS radius I and the scale 
factor a, 

L mJn = 27r^-wl.6- ) (93) 
jia a 

where j\ is the first root of J\ (cf. ( 1241 and \25V). The oscillation length is inversely proportional to a. 
Graviton oscillations become observable when the oscillation length is at least of the order of the arm 
length of a GW detector. For the ground-based interferometric detectors this requirement corresponds 
to L ~ 10 3 m. Combining this with the constraint on the maximal size of the AdS radius I from sub- 
millimeter tests of gravity, I < 10~ 4 m |42|, we find an upper limit on the warp factor a < 1CP 7 for the 
oscillation length to be detectable. Inserting this into the ratio of the amplitudes d89l and d90i we find 

A-/A+ - a 2 < 10~ 14 . (94) 

Therefore, the amplitude of a wave originating from a source on the ("hidden") £_-brane with oscillations 
which are sufficiently long to be detectable, is strongly suppressed by a damping factor a 2 as compared to 
a GW stemming from a source on the E + -brane itself. A strongly oscillating wave has to be generated by 
a source of 14 orders of magnitude stronger than that of a weakly oscillating one in order to be of the same 
magnitude, which at first sight makes the detection of RIGO's impossible. We will, however, see in the 
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next subsection that the damping factor could easily be overcome if the RS model is considered as a model 
for a putative M-theory realization of large extra dimensions. 

Beforehand we want to summarize that the main differences between the graviton oscillations consid- 
ered in this paper and the more traditional neutrino or kaon oscillations: 

First, as is readily inferred from Eq. d52l > or J54> and d55l >. respectively, the transformation between the 
gravitational modes h + and h~ living on certain branes and the graviton modes ho and hi possessing fixed 
masses is not an orthogonal (and not a unitary) transformation. This is connected with the fact that — 
in contrast to the mixing for neutrinos or kaons, where only the mass matrix is non-diagonal if written in 
terms of observable states — in the case of gravity in a two-brane world the whole kernel of the action is 
non-diagonal and, moreover, is diagonalizable only on-shell (cf. Sec. 4). 

The second and, perhaps, more important distinguishing feature of RIGO's is the dependence of the 
mixing parameters and the mass of the massive mode on the parameter a, i. e. their dependence on the 
radion. If one considers a model in which the distance between the branes is time-dependent, the mixing 
parameters and all the oscillatory effects become functions of time depending on the particular features of 
the model under consideration. 

6.4 High-amplitude RIGO's on the S + -brane from M-theory 

The bad prospects for the detection of RIGO's found at the end of the last section may considerably be 
improved if one considers the RS two-brane model as a toy model for a yet-to-be-constructed model of 
the strong-coupling limit of heterotic M-theory with one large extra dimension. Although such a model 
is lacking, we can extrapolate generic features of M-theory compactifications with small extra dimensions 
and draw some reasonable conclusions for the hidden sector of the model. These conclusions suggest the 
possibility that the hidden brane could in fact be the source of GW's with very large amplitudes. 

Cosmological M-theory models consist of an 1 1 -dimensional spacetime which is compactified on an 
S 1 /I2 orbifold with a stack of D-branes on each orbifold-fixed plane 1 431 1441 IT71 . Further 6 of the 
10 dimensions are compactified on a Calabi-Yau three-fold, while the other dimensions parallel to the 
stacks of branes, including the timelike direction, remain uncompactified. The compactification scale of 
the Calabi-Yau manifold is generally assumed to be smaller by one or two orders of magnitude than the 
size of the compact 11th dimension. Each stack of branes hosts one super-multiplet of E% gauge fields. 
This construction can become a suitable model for our universe if one of the set of Eg fields is broken 
into the subgroup Eq and further into SU(3) x SU(2) x U(l). These fields are then called the visible 
sector of the model and the corresponding stack of branes is identified as our universe, while the other Eg 
is called the hidden sector. In order to obtain a viable phenomenology it is necessary to break the Eg of the 
visible sector by gaugino condensation in the hidden sector and then mediate the influence of the gaugino 
condensate to the visible sector by its coupling to moduli fields on the Calabi-Yau manifold. 6 The volume 
of the Calabi-Yau space V at the visible sector is assumed to be of the size of the standard-model grand- 
unification scale V- 1 ' 6 ~ 10 16 GeV. In generic M-theory compactifications, it has been found necessary 
in this context to choose the volume of the Calabi-Yau space at the hidden sector smaller than at the visible 
sector, which will also lead to a stronger gauge coupling in the hidden sector than in the visible one 1441 
(cf. also |21 1). Also the volume of four-dimensional slices parallel to the uncompactified dimensions of 
the stacks of branes will decrease towards the hidden sector. It has been observed that in order to obtain a 
viable phenomenology one has to assume that the Eg on the hidden sector is also broken to a smaller group 
146) . 

In order to make connection of these M-theory setups to the Randall-Sundrum two-brane model, one 
identifies the visible stack of branes with the £ + -brane and the hidden stack of branes with the £_ -brane: 
the orbifolded coordinate of M-theory is identified with orbifolded bulk coordinate of the RS model, and 

5 Cf. ED for a non-technical description of the cosmological aspects of this class of models. 

6 Cf. [45], Chap. 18.3, p. 366-371, for a pedagogical description of this mechanism in the weak-coupling limit of Eg X Eg 
heterotic string theory. 
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the Calabi-Yau space is neglected in the effective five-dimensional description due to its small volume 
l2l . The major deviation of the RS model from M-theory ideas is that the size of the orbifolded bulk 
may be large. In M-theory models the scale factor for the uncompactified coordinates corresponding to 
the warp factor a(y) in the RS model is tied to the Calabi-Yau volume. A decreasing warp factor also 
leads to a decrease of the Calabi-Yau volume. In the M-theory solution of 1 17 1 the relation between the 
four-dimensional scale factor and the Calabi-Yau volume reads, for example, 

V(y)~a 6 (y), (95) 

where y denotes the coordinate of the orbifolded dimension. In scenarios with a large orbifolded dimen- 
sion, the dependence of the Calabi-Yau volume on the scale factor should be weaker because the string 
scale sets a lower limit for the size of the Calabi-Yau manifold. Nevertheless, if M-theory can be extrap- 
olated to setups with a large bulk, one has to expect a considerable decrease of the Calabi-Yau volume 
from the visible to the hidden brane. The crucial observation for our argument is now that the vacuum 
expectation value of the gaugino condensate rj located on the hidden brane depends strongly on the inverse 
size of the Calabi-Yau space 1471. 



-Wtt^-(S - /3T) 
boa 



(96) 



Here, a is the 10-dimensional gauge-coupling constant, bo is the beta-coefficient of the gauge-coupling 
renormalization group, S and T are moduli of the Calabi-Yau space and (3 describes loop-corrections to 
the moduli fields. From the structure of j96t we find that even a moderate decrease of V with the warp 
factor leads to a strong increase of the vacuum expectation value of the gaugino condensate, which will 
therefore be very large in a RS-like model deduced from M-theory. 

Turning now to the theory of cosmic strings, it is well known that the quadrupole moment of a cosmic 
string is proportional to the square of the vacuum expectation value rj of the corresponding symmetry 
breaking, 

/1 ~ T V 2 /uj 3 , (97) 

where T w 50 ... 100 is a numerical coefficient depending on the trajectory and shape of the string loop and 
us is the characteristic frequency of string oscillations |48 1. Therefore a large gaugino vacuum expectation 
value rj would lead to a very large quadrupole moment for cosmic strings that are produced during the 
breaking of the hidden Eg, which in turn will generate high-amplitude gravitational waves on the hidden 
brane. The described effect may easily compensate the damping factor J94I and lead to strong gravitational 
waves on the hidden brane and, as a consequence, also on the E + -brane. These in turn can be distinguished 
from waves originating from the S + -brane by detectable RIGO's. 

Note that the described mechanism is not spoiled by the dependence of the effective gravitational cou- 
pling constant on the volume of the Calabi-Yau space, G5 ~ Gu/V, because the dependence of the 
gaugino vacuum expectation value on the size of the Calabi-Yau space is much stronger than that of the 
effective five-dimensional gravitational coupling. 

At first glance our result about a strongly enhanced symmetry-breaking vacuum expectation value on 
the hidden brane compared to the visible brane seems to be in sharp contrast to the conventional interpre- 
tation of the solution to the hierarchy problem in the RS-model (cf. e.g. G||49][23|). This conventional 
interpretation states that the vacuum expectation value of a symmetry breaking on the hidden brane should 
be suppressed by two powers of the warp factor compared to the same value on the visible brane, 

V- = aV • (98) 



However, already Randall and Sundrum pointed out that the generation of the hierarchy could equally 
well be described by considering the gravitational scale on the E_-brane as the derived scale, whereas 
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the symmetry breaking scale on both branes is to be taken the same |2|. That this interpretation is in 
fact the correct one has first been demonstrated in 1 50 1 by considering the fall-off properties of on-brane 
Green functions in the Euclidean domain and has been confirmed by an analysis of the measuring process 
for masses on the E_-brane in |51 1. Therefore, all effects of the hierarchy generation in the RS model are 
already contained in the damping factor (I94> and there is no geometrical suppression of symmetry-breaking 
scales on the hidden brane in the RS model. 



6.5 Graviton oscillations on the £_ -brane 

Graviton oscillations on the £_ -brane can be treated along the same lines as RIGO's on the S + -brane have 
been treated in Secs. l6.ll and l6.3l Abandoning once again the matrix notation we find for h~ with the help 
of the spectral representation J30ft the equations of motion 
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16ttG 4 



□ - mj 



a 2 , 



a 2 T~ 



(99) 



Considering the same set of sources ( I74i we now study the gravitational waves seen by an observer on the 
negative tension brane, 



h-\T+] = Aa 2 e 



h~[T-] = Aa 2 e 



E 



i 



f Jzllrrii/a] 

h ) 



(100) 



(101) 



For simplicity we suppose again that only the massless and the first massive modes are generated, i. e. 
mi < uj < m.2- Then, for the waves generated by the source on the S + -brane and the S_-brane we have, 
respectively, 



h-[T+] = Aa 2 e- luJt (e iulr - je'v / " ! -™? r J 
h~[T-] = Aa 2 e' tut (a 2 e^ r + e V" 2 -™?^ , 



(102) 
(103) 



where A = 4 G4 fj,/r is the amplitude of the massless mode generated on the S + -brane. The amplitudes 
detected by a gravitational-wave interferometer are given by the absolute values of dl02t and dl03> . 



\hr[T+]\ =A^ 



\h-\T-]\ =A~ 



1 - 



(J2 - I) 2 

4a 2 
(1 + a 2 ) 2 



2 . 2 
sin 



sm 



/7rr \ 



1/2 
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(104) 



(105) 



Here L, the oscillation length of the amplitude modulation of the gravitational wave, is given by iU. The 
pre-factors of the amplitudes (I104i and dl05> are 



A + = I 1 ) ^A pb 1.5 of 1 A , 



A- 



J2 

(l + a 2 )a 2 Att a 2 A. 



Therefore, for the ratio of the typical intensities measured on the £_ -brane we have, 



(106) 
(107) 

(108) 
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For waves on the £_ -brane we thus find that the intensities of waves produced on this brane and on the 
other brane are of the same order of magnitude. Therefore at first glance the situation of an observer on the 
S_-brane seems to be favorable for observing graviton oscillations. 

Unfortunately, this effect is spoiled by the dependence of the graviton mass — and therefore of the 
oscillation length — on the position of the brane. In order to find the oscillation length actually measured 
by an observer on the brane, we have to use the coordinate system in which an observer does his measure- 
ments. The coordinates of an observer will always be chosen with respect to the Lorentz metric rj^ (cf. 
1511 1. No rescaling is necessary for studying effects on the S + -brane because there the induced metric 
coincides with the Minkowski metric. On the £_ -brane the "physical" coordinates used by an observer 
are related to those given with respect to the induced metric g^ v = a 2 rj^ v by 

In the equations of motion Jl 81 . the flat space d' Alembertian has to be replaced by the one with respect to 
the physical coordinates, 

To be precise, we would also have to rescale the metric perturbation h^u (with respect to the background 
metric g M „) to the measured one = h liu /a 2 {y) (with respect to the Lorentz metric). We will not per- 
form this step because we are comparing amplitudes measured on the same brane. Therefore this rescaling 
does not change the ratio between the amplitudes. In physical coordinates, the source-free equations of 
motion on the £_ -brane are given by 



phy S _ 
u phys'' — "i 



□ 



phys 



<T=o, an) 



whereas the equations of motion on the E + -brane remain unchanged. Therefore the graviton mass appears 
different when observed on a different brane, 

ml=m i =j i j, m i = — = —. (112) 

This rescaling is characteristic of bulk fields in the Randall-Sundrum two-brane model [51 50 1 ? 

Whereas the measured graviton masses mf on the S + -brane depend on the warp factor, the masses 
ral on the £_ -brane are independent of the inter-brane distance. From these masses we can derive the 
oscillation length observed on the branes for a GW with a wave-number u> > m~ . In this limit we have 

L=— . (113) 

On the £_ -brane this becomes 

L- = — = 2ir-. (114) 

m x ji 

By means of this expression, the upper experimental limit on the AdS radius I < 10~ 4 m also yields the 
upper limit on the oscillation length on the £_ -brane. However, an oscillation length of sub-millimeter 
size is clearly unobservable. Therefore graviton oscillations will be hidden to an observer on the £_-brane. 

To conclude, the situation is unfavorable for the detection of RIGO's by an observer on the £_ -brane: 
albeit being connected with strong GW's, the observer cannot detect the oscillations because their oscilla- 
tion length is far too short to exhibit itself in interferometric detectors. 



7 Note that the rescaling does not apply to on-brane fields and therefore does not affect the symmetry breaking scale considered 

■ Sec.KTSIrcf. I5T1I5BII 
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7 RIGO's in bi-gravity models 

Bi-gravity models are braneworld setups in which the effective four-dimensional gravity on the brane which 
is considered as our Universe is mediated by one massless and one light massive gravitational mode which 
have approximately equally strong couplings to matter on the visible brane. The rest of the KK spectrum is 
separated by a mass gap and its coupling to matter on the visible brane is strongly suppressed. Meanwhile 
there are various realizations of these setups |5]|9), although only the six-dimensional ones do not seem to 
suffer from a VvDZ discontinuity for the massive mode or from a phenomenologically unrealistic negative 
curvature of the visible brane. For our consideration we do not need to consider a specific realization of 
these models but can content ourselves with the general structure of their equations of motion and the 
assumption that phenomenological constraints from the VvDZ discontinuity are somehow circumvented. 
Then the four-dimensional effective equations of motion for the transverse-traceless sector of a generic 
bi-gravity model reads 



h™ = -16ttG 



□ 



Cl 
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□ ret - mi 
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(115) 



where 



Co PS ci , 
mi <C Am„ 



c„ < c , 

Am n = m n - m„_i 



n > 2, 
n > 2 



(116) 
(117) 



we have only included sources T V1S on the visible brane. The effective gravitational constant G does not 
coincide with the usual effective four-dimensional gravitational constant. Rather we have 



G 4 m G(c + ci) 



(118) 



From l^Kjl we infer that we have a strong mixing between the massless mode and the first KK mode. 
This would make bi-gravity models a natural candidate for strong RIGO's. However, this expectation is 
spoiled by phenomenological constraints on this class of models. The major experimental constraint is 
established by precision measurements of Newton's law on the orbital motion of the planets. Our treat- 
ment will closely follow that of 1 53 1 where a purely massive gravitational interaction is considered. It is 
sufficient to consider only the influence of the first KK mode, since the heavier modes will only contribute 
to the short-range dynamics. The non-relativistic static gravitational potential V(r) = —G^M/r then gets 
modified by a Yukawa potential to the form 



V(r) 



GM 



c + ci exp — 



A, 



(119) 



where M is the mass of the central body. The gravitational acceleration of a test body in the modified 
potential is g = — e r /i(r)/V 2 where e r is the unit-vector pointing from the central body to the test body 
and 



At(r) = GM [co + ci (1 + r/Ai) exp (-r/Ai)] 
GM co + ci-^ [ ^- ) +0 



(120) 



For pure Newtonian gravity fx — G4M = constant. Its value can be determined to a high precision from 
the orbit of the earth around the sun. If a massive gravitational mode contributes, the /i determined for 
the orbits of other planets will differ from the value of the earth orbit. For a planet with a semi-major 
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axis a p and a period T p , Kepler's third law yields /x 1 / 3 (a p ) = a p ^ 3 '(2tt /T p ). Therefore it is convenient to 
parametrize the deviation of a planet's fi(a p ) from the one of the earth fJ,(a®) through a small parameter 

_ ( n(a p ) x 



1 + ^= ^ . (121) 

Using Oil in dl2H . we obtain a lower bound on the wavelength of the first KK mode in terms of the 
experimental limits on 

Al= ( ci[°e- a p] V /2 _ (122) 
\^6?7p(co + ci)J 

Restricting to the case with cq m c\ w 1/2, we find 

1/2 

Ax = I I . (123) 




The most restrictive bound on rj p comes from the measurement of the Mars orbit for which we have 
rj m < —6.5 x 1CP 10 |55|. This leads to the requirement 

mi = 2tt/Ai < tt/(10 12 km) (124) 

for a realistic bi-gravity model. One might try to infer more stringent upper bounds on the mass of the first 
KK mode by considering the motions of galaxies and clusters of galaxies. However, these considerations 
rely heavily on the assumed amount of dark matter in the universe. In view of the fact that the cosmology 
of braneworld models is still in its infancy and can provide us with quite surprising models of dark matter 
(see e. g. |52|), we resist from using such constraints which would lead us to an upper mass for the first 
KK mode close to the bounds on the graviton mass given by the Particle Data Group, m g < (10 19 km) -1 

eh 

However, even the less restrictive solar-system limit J124I on mi renders RIGO's in bi-gravity models 
unobservable. Inserting J124I into the formula for the oscillation length i92\ . we find even for GW's of a 
frequency of 10~ 4 Hz, the lowest frequency observable by LISA, an oscillation length of more than 300 
light years. An amplitude modulation of that length will clearly remain undetectable. The limit on the 
wavelength of the first KK mode does however still allow for the possibility of detecting ultralight KK 
modes by studying the propagation speed of gravitational waves with LISA (cf. 1531 ). 

In addition to the bounds on the length of RIGO's from solar-system dynamics, it is interesting to note 
that some models even predict a minimal oscillation length beyond the present-day Hubble radius of our 
Universe. This is in particular the case in the Kogan-Mouslopoulos-Papazoglou version 1 9 1 of the Karch- 
Randall model |7 1 with two AdS4 branes embedded in AdSs. This model is particularly attractive because 
the longitudinal polarization of the KK gravitons is supressed in the propagator by the geometry of the 
setting. Therefore, the VvDZ discontinuity is absent I9l ll2l . 

In this model, the strongest constraint on the length of RIGO's comes from the requirement that the 
observable cosmological constant of our Universe, i. e. the effective cosmological constant A4 of one of 
the AdS4 branes, should be compatible with the observed magnitude of the cosmological constant, 

|A 4 | < nn 120 ii4, (125) 

where Mpi denotes the four-dimensional Planck mass. This requirement on the geometry of the model is 
reflected in the coefficients of the harmonics of the bulk spacetime and thereby restricts the mass of the 
first KK-mode m 1 to a fraction of the inverse size of the Hubble radius of our Universe, 

mi ~ e" 135 /^ 1 ) . (126) 
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Using the expression for the lower limit of the oscillation lenght ( 19 U . the value for the first KK mass leads 
to an oscillation length for RIGO's, 

imin ~ e 135 ^ 1 » H- 1 , (127) 

which is much larger than the size of the horizon of our Universe and thus renders RIGO's unobservable 
in principle in this particular model. 

8 Conclusions 

We have studied the phenomenon of radion induced graviton oscillations in the two-brane world. First, we 
have established a method to extract the particle content from the nonlocal braneworld effective action of 
1 26 1 for the Randall-Sundrum two-brane model. This method is easily extendable to arbitrary two-brane 
models and provides the missing link between nonlocal braneworld actions, which are specifically suited 
for the treatment of cosmological problems 1 26 , 56 , 57 1, and spectral representations of the effective action 
particularly suited to particle-phenomenology considerations 1581 1511 . 

From the equations of motion for the graviton sector, we have found a mixing of massless and massive 
modes which depends parametrically on the radion. This mixing leads to the effect of radion induced 
graviton oscillations which is in principle observable. RIGO's are a feature of every higher-dimensional 
spacetime model, since there will always occur amplitude modulations in GW's which are a mixture of 
a massless mode and KK modes. However, in traditional models with flat extra dimensions, the mass of 
the first KK mode is so big that it will neither be produced by astrophysical sources nor lead to oscillation 
lengths of macroscopic size. In contrast to this, warped geometries allow KK mode masses which are so 
low that they can lead to oscillations of detectable length. In particular, waves from sources on the hidden 
brane show strong oscillations on the visible brane. The amplitudes of these waves are strongly suppressed 
by the warped geometry and thus, at first sight, seem to remain undetectable. However, by using simple 
M-theory motivated scaling arguments, we have demonstrated that one should expect a network of cosmic 
strings on the hidden brane which would produce a background of high amplitude gravitational waves. The 
M-theory scaling properties considered may easily compensate the geometrical damping of GW's from the 
hidden brane. Therefore GW's with RIGO's stemming from the hidden brane may actually become a 
relevant effect for gravitational-wave astronomy. The characteristic pattern of RIGO's may even help to 
discriminate between noise and the stochastic GW background in GW interferometers. 

At first sight another natural candidate for observable RIGO's are bi-gravity models. In these models the 
massless graviton and the lightest massive one are coupled to matter with nearly equal strength and, thus, 
produce strong oscillations. Unfortunately, bi-gravity models which exhibit RIGO's of experimentally 
detectable oscillation-length are already ruled out by precision measurements of Newton's law on solar- 
system scales. 

The effects connected with the mixing and oscillations of quantum states in a multidimensional space- 
time have received some attention in the literature [59 60 , 61 62 1. Neutrino mixing and oscillations were 
reconsidered in l59*l l60ll6Tl . while the mixing between quarks and an attempt to explain the origin of CP 
violation was described in 1621 . Nevertheless, to our knowledge, the possibility of oscillations between 
different graviton states and their potential observability in GW interferometers has not been considered 
before. 
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A Diagonalization of the kinetic and mass terms 



Here we describe in detail the diagonalization of the kinetic and massive operators and the calculation of 
the mixing parameters from Sec. [5] The procedure of the diagonalization of the matrices D and M will 
include three successive operations. First we shall find the eigenvalues Ai, A2 of the operator D and shall 
construct an orthogonal matrix O which provides the transition from an old basis to the basis of normalized 
eigenvectors of the operator D. Application of the transformation O to the matrix D transforms it to a 
diagonal matrix whose elements coincide with the eigenvalues of D, 



ODO T = 



Ai 
A 2 



(128) 



Secondly, in order to transform this matrix into the unit matrix, we act on the eigenvectors by the 
generalized dilatation matrix A, 



A = 



AT 








Of course, the action of the inverse operators A 1 on the matrix ( 11281 yields the unit matrix, 
A- 1 OT>0 T A- 1 = 



The simultaneous action of the operators O and A on the matrix M J62I transforms it into 

M = A- 1 OMO T A-\ 



(129) 



(130) 



(131) 



which is still degenerate. The non-zero eigenvalue of this matrix, which gives the squared mass of the 
massive graviton, is equal to the trace of this matrix, 



M = TrM. 



(132) 



As a third step, in order to finally obtain an expression for the transformation from the initial gravita- 
tional modes h + and h~ to the new modes Km and ho corresponding to the massive and massless gravitons, 
we should use another orthogonal rotation Q, diagonalizing the matrix M in such a way that its diagonal 
elements coincide with its eigenvalues, i. e. 



QMQ J 



M 2 




(133) 



Of course, the application of the orthogonal transformation Q to the transformed kinetic matrix dl30> does 
not change it because it is a unit matrix. 

Thus, one can formulate the transformation between the modes h + and h~ and the eigenmodes of the 
operator (16 1 1 in the form 



h M 
h 



QAO 



h+ 
hr 



(134) 
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The inverse transformation is 

O t A- 1 Q t 



h+ 
h~ 



h M 
h 



(135) 



The transformations (II 34i and (I135> are not orthogonal. This is connected with the fact that we had to 
apply the generalized dilatation matrix A to normalize properly the kinetic part of the effective action. 
The eigenvalue equation for the kinetic matrix D, defined in Eq. d63l >. has the following form, 



A 2 - Aa(a 4 + 6a 2 + 1) + 3a 2 a 2 (a 2 + l) 2 = , 
where we have introduced the abbreviation a = (1 — a 2 )/3a 2 (l + a 2 ) 2 . The solutions are 



a(a 4 + 6a 2 + 1) oVa 8 + 14a 4 + 1 

Ai 9 — i . 

2 2 
Now, we substitute these eigenvalues into an eigenvector equation, 

(D - A 1>2 I)^i, 2 = 0, 



(136) 



(137) 



(138) 



where I is a unit matrix and ip± t 2 are eigenvectors corresponding to the eigenvalues Ai,2, respectively. We 
choose the normalized eigenvectors 



i 



Dl2 

Xi-Du 



(139) 
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D 12 
A 2 - D n 



(140) 



where Dy are the corresponding elements of the matrix D, while the normalization factors are equal to 



Pn-Ai) 2 , 



N2 = \/Dl 2 + (D 11 -\ 2 )< 



(141) 



Correspondingly, the orthogonal matrix providing the transformation from the old basis 
the new basis ifii , ip 2 (given by the expressions ( I139> , ( 11401 ). has the form 

D 12 D 12 



' 1 " 




' " 
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to 
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Ni N 2 
Ai — Du X 2 — Dn 



Ni 



N 2 



(142) 



Using the expressions for the orthogonal matrix O and for the generalized dilatation matrix A, one can 
get an expression for the rotated matrix M by substituting Eqs. dl42> and ( 11291 into Eq. 11311 . 



M = A' 1 OMO T A^ 1 = (3 



M 2 
Mm 



Mm 



^/XlX 2 

~ 2 
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VA1A2 A2 



(143) 
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where 



M = a 2 D 12 f — ) , (144) 

N\ N 2 



and 



a 2 (r — s) r + s 



Ai — A2 



No 



D u - D 2 



(145) 



(146) 



Now we are in a position to calculate the value of the mass of the massive graviton mode by substituting 
the formulas dl43> - dl46> . together with formula dl 371 and the explicit values of the matrix elements 
Dij from formula ( 16 H . into Eq. (1 1 32b . Straightforward but rather cumbersome calculations lead us to the 
simple expression 



M 



24a 2 (1 



Z 2 (l-a 2 ) 2 



(147) 



To accomplish the task of simultaneous diagonalization of the matrices D and M, we should find the 
matrix Q rotating the matrix M to the diagonal form dl 33l >. One can find normalized eigenvectors of the 
matrix M by solving the equations 



(M - M 2 I)</>! = , 
M0 2 = 0. 

These eigenvectors have the form 



(148) 
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where 



N = \J\ 2 M 2 + Aim 2 . 
Correspondingly, the orthogonal matrix Q reads 
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N 



X 2 M 
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v/A^M 



(149) 



(150) 



(151) 



Substituting the matrices O, A and Q, given by Eqs. (I142L J129l > and j!51> . into Eq. dl34l i. we find the 
expressions ( I67i . i68l , describing the transformation from the old graviton modes h + and hr to the modes 
Km and ho, which have already been given in Sec. [5] Analogously, one can also obtain the expressions 
(j69jl, J70> presented in Sec.|5]and describing the inverse transition from the eigenmodes of the diagonalized 
Hamiltonian hu and ho to the old graviton modes h + and h~ . 
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